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Abstract 

We report our numerical lattice QCD calculations of the isovector nucleon form factors for the 
vector and axialvector currents: the vector, induced tensor, axialvector, and induced pseudoscalar 
form factors. The calculation is carried out with the gauge configurations generated with Nf = 2+1 
dynamical domain wall fermions and Iwasaki gauge actions at (3 = 2.13, corresponding to a cutoff 
a -1 = 1.73 GeV, and a spatial volume of (2.7 fm) 3 . The up and down quark masses are varied 
so the pion mass lies between 0.33 and 0.67 GeV while the strange quark mass is about 12 % 
heavier than the physical one. We calculate the form factors in the range of momentum transfers, 
0.2 < q 2 < 0.75 GeV 2 . The vector and induced tensor form factors are well described by the 
conventional dipole forms and result in significant underestimation of the Dirac and Pauli mean- 
squared radii and the anomalous magnetic moment compared to the respective experimental values. 
We show that the axialvector form factor is significantly affected by the finite spatial volume of 
the lattice. In particular in the axial charge, gA/dv, the finite volume effect scales with a single 
dimensionless quantity, m n L, the product of the calculated pion mass and the spatial lattice extent. 
Our results indicate that for this quantity, m n L > 6 is required to ensure that finite volume effects 
are below 1%. 

PACS numbers: 11.15.Ha, 11.30.Rd, 12.38.Aw, 12.38.-t 12.38.Gc 
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I. INTRODUCTION 



The isovector nucleoli form factors are probes for nucleoli structure associated with the 
isovector vector and axialvector currents, = u^^d and = U'-f^^d, with up- and down- 
quark spinors u and d. From these currents, four isovector form factors arise in neutron j3 
decay: the vector (F v ) and induced tensor (Ft) form factors from the vector current, 



(p\V+(x)\n) = u p hpF v {q 2 ) + a^q x F T (q 7 



u n e 



iq-x 



(1) 



where Fy is equivalent to F\ and Ft to F2/(2Mn) in the isovector part of electromagnetic 
form factors under the isospin symmetry, and the axial (Fa) and induced pseudoscalar (F P ) 
form factors from the axialvector current, 



{p\A+(x)\n) =u % 



iq-x 



(2) 



ltilhF A (q )+iq^F P (q 

We use the Euclidean metric convention as in the recent RBC works [l, 2]. Thus q 2 stands 
for Euclidean four-momentum squared, and corresponds to the time-like momentum squared 
since q 2 M = —q 2 < in Minkowski space. Here q = p n ~Pp is the momentum transfer between 
the proton (p) and neutron (n). 

The vector-current form factors have been studied experimentally with high accuracy 
at both small (< 1 GeV 2 ) and large (> 1 GeV 2 ) momentum transfers, through electron 
elastic scattering off proton and nuclei [31. Early experiments revealed that the proton 

nnnQ 

is a composite particle [4], [5J, [6J, [7[ : i. e. non-zero Pauli and Dirac mean-squared radii and 
anomalous magnetic moments were measured among other observables. Recent experiments 



have improved the accuracy of these 
QCD predictions have been observed 



brm factors and deviations from earlier perturbative 
3- 



As is well-known, the isovector axialvector current is strongly affected by the spontaneous 



chiral symmetry breaking in the strong interaction [101. A consequence for the nucleon 
is that the isovector axial charge qa deviates from the corresponding vector charge gy. 
These isovector vector and axialvector charges, respectively the vector and axialvector form 
factors at the zero momentum transfer, are most accurately measured in neutron beta decay 
experiments: g A /gv = F A (0)/F v (0) = 1.2695(29) (l2|. Whether lattice QCD calculations 
can accurately reproduce this ratio, gA/gv, is an important test of lattice QCD. 

The spontaneous breaking of chiral symmetry also means that the corresponding form 
factors of the axialvector current are strongly coupled with the Nambu-Goldstone particles, 



i.e. the pions. Using the axial Ward-Takahashi identity and the pion-pole dominance as- 



sumption on the induced pseudoscalar, one can derive the Goldberger-Treiman relation [13] , 
which relates the nucleon mass (Mjv), the axial charge (^a), the pion decay constant (F n ) 
and the pion-nucleon coupling (g^NN)'- M^9a — F-k9-kNN- It is an interesting challenge for 
lattice QCD if it can reproduce this relation. 



The q 2 dependence of the axialvector form factor has also been studied in experiments 



It again provides a stringent test of QCD through a comparison of lattice QCD calculations 
with such experiments. While recent experiments report the induced pseudoscalar form 



factors 



15 



16], it is less well known than the other form factors. Hence this provides an 



excellent opportunity for lattice QCD to play a leading role and guide future experiments. 

In the past years, many lattice QCD studies have been made for these isovector form 
factors in the above-mentioned contexts Q]- Many earlier works Q, Ji], 



20 



21 



221 ] were 



performed either in the quenched approximation, neglecting dynamical sea-quark effects 
or were either limited to two dynamical flavors of Wilson fermion quarks that explicitly 



violate chiral symmetry 



23 



24 



quarks [26, 



27 



28 



251 1 . limited to non-unitary combination of valence and sea 

n 

29j . or just two dynamical flavors of domain- wall fermions (DWFs) J2J. 



There has also been an increasing amount of interest in the form factors of other baryons 



30 



31 



32 



27|, 



33]. 



In this paper we present our results with more realistic "2+1 flavor" dynamical quarks: 
reasonably light and degenerate up and down quarks and strange quark with a realistic 



mass are all described by the DWF scheme [34], |35j, |36|, |37j, |38] that preserves the flavor and 
chiral symmetries sufficiently Earlier studies were often performed on small spatial volumes 
(~ (2fm) 3 ) which are now widely regarded to be too small to accommodate a nucleon at 



24 



39] . We use larger spatial lattice 



light quark masses that yield realistic axial charge 
volume, as large as 2.7 fm across, to better address the finite-size question. 

The rest of the paper is organized as follows: We explain our method of calculation in 
Sec. Ull In Sec. II I II we first summarize the numerical lattice QCD ensembles used for this 
work. Then we discuss in detail the known systematic errors in the relevant form factors 
calculated on these ensembles. The numerical results are presented in section IIVI Finally, 
we give the conclusions in Sec. [Vj 

Since we vary only light quark mass in our simulation while the strange quark mass is 
fixed, in the following we call the light up and down quark mass as quark mass, rrif, in the 



4 



lattice unit, unless explicitly statec 
this study were presented in Refs. [ 



otherwise. We note that some preliminary results from 
33, y, |41 1 . 



II. METHOD 



A. Two- and three-point functions 

Following earlier studies with quenched and two dynamical flavors , we define the 
two-point function of proton 



C s (t - t src , v) = \Y. e**Tr [V A (0\ X s& t)x G (0, t SIC )\0) 



(3) 



where S is the index of the smearing of the quark operator and t src is the time location of 
the source operator. The projection operator P 4 = (1 + Ja)/2 eliminates the contributions 
from the opposite-parity state for p 2 = 42|, |43j. We use the standard proton operator, 



Xs(x) = e abc {[u s a {x)] T C^d s b {x))u s c {x) 



(4) 



where C is the matrix of the charge conjugation, and a, b, c are color indices, to create 
and annihilate proton states. In order to improve the overlap with the ground state, we 
apply Gaussian smearing 44| at the source, while at the sink we employ both local and 
Gaussian-smeared operators, S = L or G. 

In this paper we measure the nucleon isovector matrix elements for the vector and axi- 
alvector currents, 



(p\V*(x)\p) = (p|u(x)7 M u(z) - d(x)^d(x)\p), 
(p\Al(x)\p) = {p\u(x)i5W(z)-d(x)j B ipd(x)\p). 



(5) 
(6) 



While we employ the local currents in most of the calculations, the point-split conserved 



38| is used for the vector charge at the lightest quark mass which will be 



vector current 
described later. 

In order to obtain the matrix elements, we define the three-point function with the current 
J and the projector V a 

1 



4 



,iq-z 



Tr 



V a (0 1 X G (X, t snk ) J„(z,t)x G (Q, tsrc) 1 0) 



A J (^) x f(t src ,t snk ,t,M N ,E(q),q) 



(7) 
(8) 



where t S nk is the sink time slice fixed as t sn k — t src = 12, and E(q) 



(p. 



The 



ellipsis denotes the higher excited state contributions, which can be ignored for long time 
separations t sn k 3> t ^> t STC . The time independent part of A J (q) is a matrix element, 
which is a linear combination of the form factors we seek. The time dependent part of 
f(t BTC , t sn k, t, Mat, E(q), q) includes the kinematical factor and the normalization of the proton 
operator which we Gaussian smear at both the source and sink. We employ the sequential 
source method to reduce statistical fluctuations, as in Ref. 20, |45|. In the three-point 
function, initial and final proton states carry q and zero momenta, respectively. This is 
because the spatial momentum should be conserved in the function as in the two-point 
function. 

The time dependence of f (t BIC , t sr ,u, t, My, E(q), q) is removed by taking an appropriate 
ratio of the three- and two-point functions 46| 



Cg (4nk — tsrc, 0) 



Cx(£snk — t, q)Cc(t — t src , 0)Ci(t sn k — tsrc, 0) 



1/2 



, (9) 



_CL(tsnk — t, 0)Ca(t — t src , g)Ci(t sn k — tsrc, Q) 

where K = M^^j2E{q){M^ + E(q)). The ratio should display a plateau from which 
the matrix element we seek is extracted. 

For each of the vector or axialvector currents, we first obtain A J (q) in eq.® which is a 
linear combination of the form factors. For convenience, using the ratio R we define 



KM 

for the vector current, and 

A£(g,f) 



M N (M N + E(q)) 
_l(R% 3 (q,t) 
2 



ftV 5 3 



(q,ty 



iq 2 M, 



N 



iqiM. 



N 



M N (M N + E(q)Y 



(10) 

(ii) 

(12) 
(13) 



for the axialvector current. Here we also define q 2 = 2M]s[(E(q) — M/v), and P53 = (1 + 
74)7573/2 implies the z-direction is chosen as the polarization direction in our calculation. 
In the plateau region of A J (q,t) we determine the matrix elements of each current, A J (g) 
which has the following relation to the form factors: 



K(q) = F^q 2 ) 
A£(?) = F 1 (q 2 ) + F 2 (q 2 



q* 
4M 2 N 



F2(q 2 ), 



(14) 
(15) 
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for the vector current, and 



A i (" = F ^ 2) ~ MJm Fp(q% (16) 

A£(,) = AfwFpfe 2 ), (17) 

for the axialvector current. In the following we use the isovector part of the Dirac and 
Pauli form factors, F 12 , rather than the vector and induced tensor form factors. They are 
identical through the isospin symmetry except the normalization of the Pauli form factor, 
F 2 = 2MnFt- We will see that the signal of these combinations is reasonable in Sec. IIVBI 
Finally respective form factors are obtained by solving the sets of linear equations, (I14j) and 
(USD, or ([ED and (ED, at fixed q 2 . 



B. Double source method 



We find the ensemble with the lightest quark mass of rrif = 0.005 is much noisier than 
the ones with heavier mass values: it is insufficient and takes enormous amount of calcula- 
tion time to obtain reasonable statistical error if we used only a single nucleon source/sink 
combination per configuration. 

Fortunately, the time extent of the lattice, 64 x a = 7.3 fm, is very large compared 
to the inverse of the nucleon mass, M^ 1 = (1.15 GeV)" 1 = 0.17 fm. Hence, we can easily 
accommodate a pair of source/sink combinations on each configuration without letting them 
interfere with each other if the sources are separated by 32 units, as shown in Fig. [TJ We 
call this the double-source method. 

The three-point functions are calculated with the sequential source method, and the 
sink operators are placed 12 time slices from their respective sources. The number of the 
measurements is effectively doubled in this calculation while the cost remains the same as 
one single source measurement. 



III. ENSEMBLES 



A. Statistics 



The RBC-UKQCD joint (2+1) -flavor dynamical DWF coarse ensembles 47 1 a re used 
for the calculations. These ensembles are generated with Iwasaki gauge action [48| at the 
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FIG. 1: Two-point function with the double source at t = (0,32) for m/ = 0.005. 
coupling (3 = 2.13 which corresponds to the lattice cutoff of a^ 1 = 1.73(3) GeV, determined 



from the Q baryon mass 



47|. 



The dynamical strange and up and down quarks are described by DWF actions with the 
fifth-dimensional extent of L s = 16 and the domain-wall height of M5 = 1.8. The strange 
quark mass is set at 0.04 in lattice units and turned out to be about 12% heavier than 
the physical strange quark, after taking into account the additive correction of the residual 
mass, m rcs = 0.003. The degenerate light quark masses in lattice units, 0.005, 0.01, 0.02 and 
0.03, correspond to pion masses of about 0.33, 0.42, 0.56 and 0.67 GeV and nucleon masses, 
1.15, 1.22, 1.39 and 1.55 GeV. 

Two lattice volumes used are 16 3 x 32 and 24 3 x 64, corresponding to linear spatial extent 
of approximately 1.8 and 2.7 fm, respectively. The smaller volume ensembles, calculated 
only with the heavier three light quark masses, are used for a finite volume study of the 
axial charge and form factors discussed in Section HVl On the 16 3 ensembles we use 3500 
trajectories separated by 5 trajectories at rrif = 0.01 and 0.02, and by 10 at 0.03. The main 
results are obtained from the larger volume ensembles with the number of the configurations 
summarized in table [H 

On the larger volume at the heavier three quark masses, we make four measurements on 
each configuration with the conventional single source method using t src = 0, 16, 32, 48, or 
8, 19, 40, 51. At the lightest mass the double-source method is used, and two measurements 
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m f 


-^conf -^sep 


TV 

1 v mcas 


m w [GeV] M N [GeV] 


0.005 


932 a 


10 


4 6 


0.3294(13) 1.154(7) 


0.01 


356 


10 


4 


0.4164(12) 1.216(7) 


0.02 


98 


20 


4 


0.5550(12) 1.381(12) 


0.03 


106 


20 


4 


0.6681(15) 1.546(12) 



TABLE I: iV con f, iV se p and iV me as denote number of gauge configurations, trajectory separation 
between each measured configuration, and number of measurements on each configuration, respec- 
tively, on (2.7 fm) 3 volume. Table also contains the pion and nucleon mass for each ensemble. 
a Total number of configurations is actually 646. We carry out extra measurements on a subset of 
these (286 configurations) to improve the statistics using different source positions. 
b Two measurements with the double-source method gives effectively four measurements. 

on each configuration are carried out using the source pairs of (0, 32) and (16, 48), or 
(8, 40) and (19, 51). We make an additional two measurements on roughly half of the 
configurations with another source pairs. This means that we make four, double-source 
measurements on almost half of the configurations, while two, double-source measurements 
are carried out on the remaining configurations. We have checked the independence of 
these measurements from each other by changing the block size in the jackknife analysis, 
e.g., treating each source/sink measurement as independent. None of these resulted in 
significantly different error estimate. Thus in the following we treat the two double-source 
measurements performed on a single configuration, one with the source pairs of (0, 32) and 
(16, 48), and the other with the source pairs of (8, 40) and (19, 51), as being independent 
of each other. 

In the following, in order to reduce possible auto-correlations at the larger volume the 
measurements are blocked into bins of 40 trajectories each, while 20 trajectories at the 
smaller volume. The statistical errors are estimated by the jackknife method. 

B. Correlation functions 

The quark propagator is calculated with an anti-periodic boundary condition in the tem- 
poral direction, and periodic boundary conditions for the spatial directions. We employ 
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FIG. 2: Measured nucleoli energies in lattice unit at my = 0.01. Estimated energies by continuum 
and lattice momenta are also plotted. 

gauge-invariant Gaussian smearing at the source with smearing parameters (N, u) = (100, 7), 



which were chosen after a series of pilot calculations, as described in Ref. |49|. For the cal- 
culation of the three-point functions, we use a time separation of 12 timeslices between 
the source and sink operators to reduce effects from excited state contributions as much as 
possible. 

To obtain the form factors at non-zero q 2 , we evaluate the two- and three-point functions, 
eqs.([3]) and ([7]), with the four lowest nonzero momenta: p = 2ti/Lx (0,0,1), (0,1,1), (1,1,1), 
and (0,0,2), corresponding to a q 2 range from about 0.2 to 0.75 GeV 2 on the large volume, 
while on the small volume we use only the smallest two momentum transfers, corresponding 
to q 2 ~ 0.4 and 0.8 GeV 2 . All possible permutations of the momentum including the positive 
and negative directions are taken into account. 

There are several choices for the definition of the momentum in the lattice calculation, 
e.g., Pi = 2n/L ■ rii, sm(2na/L ■ rii)/a, or one determined from the measured energy in the 
two-point function. Figure [2] shows that the three energies with the different momentum 
definitions reasonably agree with each other. In the following we choose the continuum 
momentum definition pi = 2n/L ■ rii, since this simple definition gives smaller statistical 
error for the energy than the measured one. 
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C. Systematic errors 



There are two important sources of systematic error: finite spatial size of the lattice and 
excited state contamination. Chiral-perturbation-theory-inspired analysis of the former for 
meson observables suggests the dimensionless product, m n L, of the calculated pion mass 
and lattice linear spatial extent L, should be set greater than 4 to ensure that the finite- 
volume correction is negligible below one percent, and the available lattice calculations 
seem to support this. While our present parameters satisfy this condition, it should be 
emphasized that such a practical criterion is not known sufficiently for baryon observables. 
It is important to check this through the present calculations, and it is indeed an important 
purpose of this work. 

On the other hand, one should adjust the time separation between the nucleon source 
and sink appropriately so the resultant nucleon observables are free of contamination from 
excited states. The separation has to be made longer as the quark masses decrease. In our 
previous study with two dynamical flavors of DWF quarks 2| with a similar lattice cutoff 
of about 1.7 GeV, we saw systematic differences between observables calculated with the 
shorter time separation of 10, or about 1.16 fm, and longer 12, or 1.39 fm: the differences 
amount to about 20 %, or two standard deviations. This would suggest that at the shorter 
time separation of about 1.2 fm, the excited-state contamination has not decayed sufficiently 
to guarantee correct calculations for the ground-state observables 41|. There is, however, a 
price to pay for the larger time separation as the nucleon correlation function suffers from 
large statistical noise at large times, especially with light quark masses. Since the hadron 
masses are much lighter in the present work than we considered previously (the lightest pion 
mass is 0.33 GeV and nucleon 1.15 GeV) we decided to use the separation of 12 lattice units, 
or about 1.4 fm. 

While it is desirable to use a longer separation, it cannot be made too long in practice 
without losing control of statistical errors. In Fig. [3] we present the nucleon effective mass 
at the lightest quark mass, rrif = 0.005. The nucleon signal begins to decay at t = 12, 
or about 1.4 fm: this is about longest distance we can choose without losing the signal, 
and hence about as free of excited-state contamination as we can achieve with the present 
statistics. As will be shown in detail in this paper, the bare three-point function signals for 
the form factors for this source-sink separation of t — 12 are acceptable. Whether this is 
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FIG. 3: Effective mass of nucleoli correlator with Gaussian smearing applied at both source and 
sink, for quark mass m/ = 0.005. 

sufficiently long separation between the source and sink to guarantee correct calculations of 
ground-state observables remains a future problem. 

IV. RESULTS 



A. Vector and axial charges 



Much of the results and discussion in this subsection have appeared in Ref . 39(] . We repeat 
them here for convenience, and to lay some of the ground work necessary for discussion of 
the form-factor results that follow. 

At zero momentum transfer the time component of the vector form factor gives the vector 
charge, gy = Fi(0). For our calculations at the heaviest three quark masses, we use the 
4-dimensional local current. As a result, the value of g^, measured from the bare -Fi(O), 
deviates from unity, and gives the inverse of the renormalization, Zy, for the local current. 
At the lightest quark mass, m/ = 0.005, we evaluate the vector charge using the point split 



conserved vector current [38|, V4 as well. This is to alleviate a problem that arises from the 
double source method described in Sec. Ill Bl Conventionally the vector charge is calculated 
from the ratio of the three-point function with the local vector current to the two-point func- 



12 



1.42 



1.4 - 



m f =0.005 



1.38 - 



1.42 



1.4 - 



m =0.01 



1.38 - 




# 



ft 







2 



4 



6 
t 



8 



10 



12 



FIG. 4: Plateaus of for m f = 0.005 (top) and 0.01 (bottom). Statistical errors are of 
comparable sizes for the two mj values despite difference in the methods. Solid lines denote fit 
results with one standard deviation. 

tion with zero momentum, as in eq.Q; a strong correlation between the denominator and 
numerator suppresses the statistical error associated with such calculations. This correlation 
is lost in the double-source calculation and results in larger statistical errors. Fortunately, 
the three-point functions of the local and conserved currents are highly correlated, even 
in this method. Therefore we evaluate the vector charge from the ratio of the three-point 
functions ffi = C%(Q,t)/C%(6,t) at m f = 0.005. Figure 1 shows that the error in this 
ratio is as small as that coming from the single source calculation at rrif = 0.01. 

A linear extrapolation to the chiral limit yields an accurate estimate of g 1 ^ = 1.3929(17), 
as shown in Fig. [5j This corresponds to Zy = 0.7179(9) in the chiral limit and agrees 
well with an independent calculation in the meson sector j^, Za = 0.7161(1), up to the 
discretization error. 

The axial charge is calculated from the ratio of the vector and axialvector form factors 
gA = Fa(0)/Fi(0). This ratio gives the renormalized axial charge since the vector and axial 
currents, and A^, share a common renormalization thanks to the good chiral symmetry 
properties of DWF, up to small discretization error of 0(a 2 ). 

The plateaus of gA computed on volume V = (2.7 fm) 3 are shown in Fig. El We checked 
that consistent results are obtained by either fitting or averaging over appropriate time slices, 
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FIG. 5: <7y l and 1/Za obtained from the pion-to- vacuum matrix element of the conserved axi- 
alvector current \^\ - 
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FIG. 6: Plateaus of g A . V = (2.7 fm) 3 and mj = 0.005, 0.01, 0.02, and 0.03, from top to bottom. 



t = 4-8, and also by fitting the data symmetrized about t — 6. The data can be symmetrized 
because the source and sink operators are identical in the limit of large statistics. We note 
that the statistics at our lightest mass is the largest we know of for comparable simulation 
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in 



f 



0.005 



0.01 



0.02 



0.03 



(2.7 fm) 3 1.073(39) 1.186(36) 1.173(36) 1.197(30) 
(1.8 fm) 3 N/A 1.066(72) 1.115(58) 1.149(32) 

TABLE II: Summary of axial charge, g^, for both volumes 
l 1 1 
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FIG. 7: Axial charge qa together with two-flavor 



action 



26 



2| and quenched 



29J] calculations. Recent N f = 2 + 1 DWF by LHP [291] is also plotted. 



DWF, and mixed 



parameters in the literatures. Results obtained from the fit using the unsymmetrized data, 
presented in the figure with one standard deviation, are employed in the analysis. These 
results are compiled in table HT1 

Figure [7| shows that the (2.7 fm) 3 data are almost independent of the pion mass (squared) 
except for the lightest point which is about 9% smaller than the others. A set of the results 
obtained with a smaller volume, (1.8 fm) 3 shows a similar downward behavior, albeit with 

Q 

relatively larger statistical uncertainties. An earlier two flavor calculation by RBC [2J with 
spatial volume (1.9 fm) 3 and 1/a— 1.7 GeV showed a clear downward behavior, but it sets 
in at heavier pion mass. 

We suspect that this pion mass dependence driving gA away from the experimental value 
is caused by the finite volume of our calculation. Similar behavior was observed in quenched 
[l|, 2(| and was predicted in a model calculation 50]. However, for pion masses 



DWF studies 
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close to our lightest point such a sizable shift is not observed when V is larger than about (2.4 



fm) 3 , not only in the quenched case, but also the 2+1 flavor, mixed action calculation in 



and their updated results 



291 ] . Both the results of quenched 



HQ 



and mixed action 



calculations on larger volumes are presented in Fig. [7J On the other hand, our results suggest 
that a volume of V = (2.7 fm) 3 is not large enough to avoid a significant finite volume effect 
on gA when m n < 0.33 GeV in dynamical fermion calculations. It is worth noting that the 
bending of the axial charge comes from only the axialvector part Fa{0), since the vector 
part i*i(0) does not have such a pion mass dependence (see Fig. [5]). 

In order to more directly compare the various results, we plot gA against the dimensionless 
quantity, m n L, in the top panel of Fig. [SJ We find that the 2+1 flavor results on both 
volumes reasonably collapse onto a single curve that monotonically increases with m n L; 
in other words, they exhibit scaling in this variable. The two flavor results display a 
similar behavior which is also evident in dynamical two flavor (improved) Wilson fermion 
calculations as shown in the middle panel 2J, |27J, |5l|] for the unitary points K sea = K, va \, with 
various volumes (0.95-2.0 fm) 3 , pion masses 0.38-1.18 GeV, and gauge couplings. While 
the trend is similar in the quenched DWF case with pion masses in the range 0.39- 

0.86 GeV and 1/a = 1.3 GeV (see bottom panel), the scaling is violated for the point with 
smallest m n L on V = (2.4 fm) 3 . The lightest point does not follow the (1.8 fm) 3 data: they 
differ by 2.5 standard deviations (a) at m n L ~ 5, suggesting that there are non-universal 
terms that depend separately on and V. In particular, this effect may be due to the 
presence of a quenched chiral log 52|. From Ref. 52[, the size of the effect at this mass 
can readily explain the discrepancy observed with the dynamical m n L scaling. Note, at this 
mass, but going to V = (3.6 fm) 3 , no finite volume effect is detected in the quenched case 
as can be seen in Fig. [7J 



The mixed action, 2+1 flavor result with a similar volume [26|, |29j, is denoted by the left 
triangle in the top panel. We plot their recent result at our lightest point 29]. At heavy pion 
masses the results are statistically consistent with our larger volume data and essentially 
independent of m n L. At m^L ~ 4.5 the mixed action result, however, is larger than ours by 
(a combined) 2.1a, and lies between our lightest result and the quenched DWF result with 
(2.4 fm) 3 volume 20j (the up triangle in the figure). 

A possible explanation of the differences is that it is simply a dynamical fermion effect 
as discussed in Ref. 39[. While the mixed action result at m^L ~ 4.5 has come down from 
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FIG. 8: m n L scaling of the axial charge. Top, middle, and bottom panels are dynamical DWF, 
dynamical Wilson, and quenched DWF results [l|, |20|], respectively. In the top panel circle, square, 



201 ] and mixed action 



26 



and diamond denote 2+1 flavor larger, smaller volumes, and 2 flavor data 2|, respectively, and 
open up and left triangles denote quenched 



291 ] data, respectively. In the 



middle panel closed symbols denot e dy namical improved- Wilson fermions 



27 



24f | . and open symbols 



5l|. In the bottom panel, the open triangle is the same as in 



are dynamical Wilson fermion 
top panel. 

higher value with larger error (the previous result was consistent with the quenched result 
at the similar m n L), the explanation using the systematic error of the partially 

quenched effect of the mixed action results might be valid in the present data. If the sea 
quark is effectively heavy, a mixed action calculation will be closer to the quenched case. 
Mixed action chiral perturbation theory reveals the presence of partially-quenched logs whose 



size is consistent with the observed effect 



55 



56(, as in the quenched theory. We should 



note that the preliminary result obtained by LHP 29(] at the same simulation parameter 
as our lightest point appears inconsistent with our result (see Fig. [7]). We will discuss the 
difference later in this section. 
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For the chiral extrapolation of gA, we attempt to include the finite volume effect in our 
data. While the pion mass dependence of gA, including the finite volume effect, has been 
investigated in the small scale expansion (SSE) scheme of heavy baryon chiral perturbation 
theory (HBChPT) [24 1. the size of the finite volume effect on V = (2.7 fm) 3 predicted in 
SSE is less than 1% in our pion mass region. The correction is much too small to account 
for the observed finite volume effect in our data. This suggests that the finite volume effect 
in HBChPT, which is estimated by replacing all loop integrals by summations, is not the 



leading finite volume effect in gA, as in the e regime [57|. We also note that our attempts 



to fit the mass dependence of the data to HBChPT failed, which is likely due to the heavier 

nnr 

quark mass points being beyond the radius of convergence of ChPT [2|, |47|, 



1581- 

Instead of the SSE formula, we assume the following simple fit form, including the finite 
volume effect in a way that respects the scaling observed in the data, 

A + Bml + Cf v (m n L), (18) 

with fy(x) = e~ x , and where A, B, and C are fit parameters. The third term corresponds 
to the observed finite volume effect, taken as a function of m n L only, and vanishes rapidly 
towards the infinite volume limit, L — > oo, at fixed pion mass. The same m^L dependence 
appears in one of the finite volume effect contributions in Ref . {^J] . We note that this simple 
form is used to estimate the finite volume effects in the data but not the value of gA in the 
chiral limit at fixed L. In the end, we choose this simplest form, in part, because the fit 
result at the physical point is not sensitive to the particular choice of fy(x), as discussed 
below. 

In Fig. [9] we see that the 2+1 flavor data are described very well by this simple fit 
(x 2 /d.o.f. = 0.57), using data computed on both volumes simultaneously. The L — > oo 
extrapolation (solid line) in turn allows an extrapolation to the physical pion mass (m n = 135 
MeV), gA = 1.19 (6) (4), where the first error is statistical. The second error is an estimate of 
the systematic error determined by comparing; this result with that from fits using different 



choices of fy(x), e.g., the full form in 59J, x 3 , and rn^e x /x 1 ^ 2 . The latter is similar to 



HBChPT when m v L » 1 U£ 



60 



6l|]. The results of some of the fit forms are summarized 
in table IHIl The extrapolated value is not sensitive to the choice of fy, and is also consistent 
with a linear fit to the three heaviest points on the larger volume, gA = 1.17(6). The present 
data are insufficient to determine the detailed form of fy, but do allow a reasonable estimate 
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FIG. 9: Chiral and infinite volume extrapolation of g& with finite volume effect fy = e x in 
eq. (|18p . Fit is carried out with data on both the volumes, simultaneously. 



fv 



A 



B 



C x 2 /d.o.f. < 



phys 



1.187(57) -0.12(14) -8.1(3.9) 0.57 1.187(55) 
1.226(70) -0.05(15) -11.3(5.2) 0.49 1.225(67) 
m 2 e -m.i / / v ^j; L i4 8 (4 6 ) -0.09(12) -75(41) 0.80 1.150(44) 



(m n L) 3 

rriTrL I 

N/A 



1.172(58) 0.05(1.7) 



0.17 1.173(55) 



TABLE III: Fit results of gA, together with the extrapolated result at m,P hys = 135 MeV. In the 
last row, the linear fit result using only the three heavier points at V = (2.7 fm) 3 is presented. 

of the finite volume effect. 

We also fit our data, with and without the fy term, to the 2-loop formula from 
HBChPT 58| and find that the extrapolated result is less than 1 and that the fits are 
generally unstable. This is due to the many unknown low energy constants which cannot 
be determined accurately from only four data points, even if some of them are fixed. More 
importantly, though the 2-loop formula extends the range of applicability of the chiral expan- 
sion, it is still only large enough to include our lightest point, as demonstrated in Ref. 58]. 
The systematic error arising from the difference of the renormalization constants for and 
Vfj, is much smaller than the quoted systematic error. From the fit result with fy(x) = e~ x , 
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we estimate that if one aims to keep finite volume effects at or below 1%, then for m n = 0.33 
GeV, spatial sizes of 3.5-4.1 fm {m n L m 5.9-6.9) are necessary. 

As mentioned, our lightest result on (2.7 fm) 3 differs from the preliminary findings from 



LHP 



29j shown in Fig. [7] by 1.8 a. These calculations are carried out with the same 



parameters except for the operator smearing and the time separation between the source 
and sink operators, At = £ sin k — t STC . So, while it is possible that this difference is simply 
due to the limited statistics in the preliminary result in 29|, there is the possibility that 
this difference is due to a systematic error stemming from contaminations of higher excited 
states. These contaminations will be negligible when the time separation of the two nucleon 
operators in the three-point function, eq. (j7jh is large enough. The large separation, however, 
causes the statistical error of the three-point function to increase. Thus, we employ a time 
separation of At = 12, as described in Sec. [Til while LHP uses At = 9. While further 
investigation of this difference is desirable, it is beyond the scope of this paper. 

Although there may be a systematic difference between our result and the result of 
LHP at the lightest quark mass on the (2.7 fm) 3 lattice, all recent results (before chiral 
extrapolation) with dynamical quarks are about 10% smaller than the experimental value. 
In order to make a precise test of (lattice) QCD with the axial charge, further study of the 
systematic errors as the quark mass is decreased towards the physical point is required on 
large volumes. 



B. Form factors of the vector current 



In this subsection we discuss the isovector part of the Dirac and Pauli form factors, F\(q 2 ) 
and F 2 (q 2 ). In Fig. [10] we present the ratios of the three- and two-point functions, and 
A^ defined in eqs. fllOl) and (II ip . at the quark mass mf = 0.01 for each momentum transfer. 
We find excellent plateaus in the middle time region between the nucleon source and sink 
operators at t — and 12 for the smaller momenta, while the plateau at q oc (2, 0, 0) is not as 
well behaved and has a larger error. This is interpreted as simply a statistical fluctuation. In 
order to remove this wiggle, we would need more statistics at this momentum. To determine 
the values of the ratios, we perform a constant fit in the time interval, t = 4-8, for all 
momentum combinations. 



20 



1.6 
1.4 



1.2 - 



1 

w 

< 0.8 



0.6 
0.4 



~\ — 1 — i — 1 — i — 1 — r 



; * i i t f f f ! !i 



0.2 « ' i- 



J i L 



J . L 



• (0,0,0) 
■ (1,0,0) 

♦ (1,1,0) 

▲ (1,1,1) 

•4 (2,0,0) 

_1 I l_ 



10 12 



14 



1.6 
1.4 
1.2 
1 



< 



0.8 h 

0.6 
0.4 
0.2 







j — • 1 • 1 • — r 



i 1 1 1 1 1 r- 



* * i " 



i 



♦ 
▲ 



(1,0,0) 
(1,1,0) 

(1,1,1) 

(2,0,0) 



i\ " it i\ 



1 i I i I i l_ 



J ± ± 



_i i i_ 



10 12 14 



FIG. 10: Ratios of 2- and 3-point functions of the vector current, and Ay, at mj = 0.01. 



The form factors are obtained by solving the linear equations (1141) and ( |T5l) . 

, for all q 

1 + r 



(19) 
(20) 



where r = g 2 /(4M^). All the values of the two form factors are shown in table [VTT1 in the 
appendix. 
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FIG. 11: The Dirac form factor, Fi(q 2 ), normalized to unity at q 2 = 0. 
1. Dirac form factor F\{q 2 ) 

Let us now turn our attention to the momentum dependence of the Dirac form factor. In 
Fig. [ED we present the form factor at each quark mass normalized by the respective values 
at zero momentum transfer. 

Phenomenologically the form factor is described by the conventional dipole form, 

F ' { " 2 > = ( i +9 W (21) 

where Mi is the dipole mass for this form factor, and fits to experimental data give M\ = 



0.857(8) GeV 12j. In order to test the dipole form using our lattice results, for convenience 



we define an effective dipole mass 



cir 



(22) 



Figure [T2l shows that the effective dipole mass at rrif = 0.01 is almost flat against q 2 . This 
means that the form factor is well explained by the dipole form eq. ( 12TI) in the q 2 region where 
we measure. The figure also shows that the effective mass is consistent with the dipole fit 
result as expected. 

We estimate the Dirac root mean-squared (rms) radius from the dipole mass obtained by 
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result 
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23j. Result of the dipole fit (solid line) with one standard deviation (dashed line) is also 



presented. 



the fit as 



2vl/2 



Mi 



(23) 



whose results are presented in table [IVl Figure [14] shows the pion mass dependence of our 
results for the rms radius. Here we also compare with other lattice calculations and the 
experimental value. Our results show a near-linear dependence in the pion mass squared 
which is quite different from the axial charge in Sec. IIV Al This suggests that the Dirac form 
factor is less sensitive to finite volume effect than g^, and this is confirmed by an analysis 
of our results obtained on a smaller volume (1.8 fm) 3 , shown in Fig. [131 The smaller volume 
results are summarized in the appendix. Our results can be fit linearly and extra jo ated 



to a value 27% smaller than experiment, 0.797(4) fm. Other lattice calculations [l|, [2|, [22], 
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25] show similar trends. The recent results of the mixed action calculation [29[ are also 



statistically consistent with our data and fit line. 

This quantity is expected to logarithmically diverge in HBChPT 



62 



63 



64) at the 



chiral limit: such a behavior will help in bringing our present extrapolated results closer 
to experiment. However, our results at m n > 0.33 GeV fail to reveal such a logarithmic 
divergence. A naive determination of the HBChPT parameters at the physical point give 
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FIG. 13: Comparison of F\ with larger and smaller volumes denoted by closed and open symbols, 
respectively, at each quark mass. 

the logarithmic contribution shown in Fig. [14] by the solid line. Future work will require 
simulations to be performed at lighter quark masses, e.g., < 0.2 GeV, if such logarithmic 
effects are to be seen in lattice results of the Dirac radius. 



2. Pauli form factor F2(q 2 ) 

Figure [15] shows the momentum-transfer dependence of our results for the Pauli form 
factor at each quark mass. These values are tabulated in table I VIII The form factor is 
renormalized by -Fi(O). 

This form factor can also be described by the conventional dipole form, 



^(0) 

(i + qyMty 



^(S 2 ) = n , T,iw ( 24 ) 



with M2 = 0.78(2) GeV and ^(0) = 3.70589 extracted from fits to experimental data. In 
contrast to the Dirac form factor, there are two parameters, the over-all strength F 2 (0) and 
the dipole mass M 2 : the former gives the isovector part of the anomalous magnetic moment, 
A*p — A*n — 1) and the latter the Pauli mean-squared radius, (r|) = 12/Mf , as in the Dirac 
case. We fit the form factor with these two parameters. 

To check reliability of the dipole fit, we measure the ratio of the Sachs electric and 
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231 ]. and mixed action 



2| and quenched DWF lj, two- flavor and quenched Wilson 
calculations, respectively. A prediction from HBChPT with the experimental result 
plotted. 
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29] 



231 ] is also 



magnetic form factors 



G M (q 2 ) K v T {q) 



(25) 



G E (q 2 ) AY(qY 

At zero momentum transfer, we obtain 1 + ^(0) from the ratio. Figure [16] shows that the 
result for Ge{q 2 )/Gm{(1 2 ) — 1 at g 2 = 0, obtained via a linear fit in q 2 , is consistent with the 
determination from a dipole fit of F 2 (q 2 ). 

In Fig. [T7]we present the anomalous magnetic moment of the nucleon, determined by the 
dipole fit presented in table IIVt together with some other lattice QCD calculations and the 



experimental value. Our present resu 
with previous lattice calculations [l[ 



ts slightly decrease with the pion mass, in agreement 
231 ] . They extrapolate well linearly in the pion mass 
squared, and result in a value 26% smaller than the experiment. This result at the physical 
pion mass is consistent with those of previous calculations [l|, [ji]] using a linear fit. 

We present in Fig. [TH] the result of the Pauli rms radius. These results are obtained from 
a dipole fit and summarized in table IIVI Some other lattice QCD calculations Q, [^J] are 
also plotted in the figure for comparison. We find the lightest point to be slightly smaller 



25 



4 



■ s 2 • m=0.005 

35 [\ F 2^ ) ■ m f =0.01 

. \ ♦ m=0.02 

\ i A m=0.03 



2.5 
2 

1.5 
1 

0.5 




T 



\ 

s 



I 



Si 

N 

s, 

V. 



- ■ experiment ~ 



i ■ ■ ■ i ■ ■ ■ i 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

q 2 [GeV 2 ] 



FIG. 15: The Pauli form factor, ^(g 2 ), renormalized by Zy = l/.Fi(0). The dashed curve is a fit 
to experimental data. 

1 I 1 1 1 I 1 1 



I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 1 1 1 I 

5h - ■;- 



4 - 



F,(q 2 ) 



G M (q 2 )/G F (q 2 )-l 



M vn ' E 

experiment 



iment j 

i T-—J' 



1 - 



o' i ■ ■ 'J, 



■ ■ ■ ■ i ■ ■ ■ i ■ ■ ■ i ■ ■ ■ i ■ ■ ■ i 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.! 

q 2 [GeV 2 ] 



FIG. 16: Dipole fit with F2(q 2 ) and linear fit with ratio of electric and magnetic form factors 
G M (q 2 )/G E (q 2 ) -latm/ = 0.01. 

than the results at the other quark masses, albeit with a large error. Thus, we consider this 
pion mass dependence is due to statistics, not a finite volume effect as in the axial charge, 
and this is confirmed by our results from the smaller volume simulations in Fig. O The 
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respectively, at each quark mass. 

results are reasonably fitted by a linear function of the pion mass squared, and we obtain 
( r i) 1//2 = 0.64(6) fm at the physical pion mass. This result again is 27% smaller than the 
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231 ] is also plotted. 



m f 



0.005 



0.01 



0.02 



0.03 



m 



phys 



(r 2 ) 1 / 2 [fm] 0.564(23) 0.548(24) 0.520(31) 0.485(16) 0.584(23) 
(rf) 1 / 2 [fm] 0.578(60) 0.690(61) 0.536(21) 0.537(38) 0.636(57) 
F 2 (0) 2.82(26) 3.40(35) 3.11(21) 3.55(19) 2.75(28) 

TABLE IV: Dirac and Pauli rms radii (r 2 ) 1 / 2 , (r^) 1 / 2 , and anomalous magnetic moment ^2(0) = 
fJ'P — — 1- The linear fit results at mP hys = 135 MeV are also presented. 

experimental value, 0.88(2) fm. 

Here again the quantity is expected to diverge as \j Jm^ in the chiral limit in 



HBChPT 
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63 



641 ]. however our results do not indicate such divergence. In contrast 



to the Dirac radius case, perhaps because of the larger statistical errors, HBChPT can si- 
multaneously fit the experiment and our data. The fit with the prediction, however, gives 
larger x 2 /d.o.f. (degrees of freedom), and twice larger value at the physical pion mass than 
the linear fit. We need further light quark mass calculation with better statistics to test the 
prediction in the lattice QCD calculation. 
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C. Form factors of the axialvector current 



In this subsection we show the form factors obtained from the axialvector currents, Fa{q 2 ) 
and Fp(q 2 ). They are extracted from the ratios of three- and two-point functions defined 
in eqs. (H2j) and (I13p . Figure I2"0l shows that the typical plateaus of the ratios with the A 3 
component of the current at m/ = 0.01 are reasonably fiat in the middle time region between 
the source and sink operators. We plot the ratios (g 3 = 0, t) and Af (g 3 ^ 0, t) separately, 
since Af(g 3 ^ 0,t) contains both form factors, while Af(g 3 = 0,t) contains only ^(g 2 ). It 
is worth noting that there is no Af (g 3 = 0, t) in the case of q oc (1, 1, 1). A^(g, t) has a slope 
in the range t = 1-8 with large statistical errors as shown in the bottom panel of Fig. [2UJ 
We consider the slope to be caused by poor statistics in the data. The values of the matrix 
elements for all the ratios are determined by constant fits with the range of t = 4-8. 

Using the relations eqs. f|T6|) and (|T7l) . the two form factors are determined through 
the following equations which depend on the spatial momentum transfer in the three-point 
function, 

A£(g 3 = 0) for n = 0,l, 2, 4 

F A{q 2 ) = { a 2 (26) 

M N (M N + E(q)) 

A$(q)/M N forn = 2,3 

W)HM^w (A . (93 = 0) _ A . (537 , 0))forn=M m 

where n = q 2 ■ (L/2n) 2 . The results for the two form factors are summarized in table IVTT1 
in the appendix. 



1. Axialvector form factor F^(g 2 ) 

Figure [21] shows the axialvector form factor at each quark mass, which is renormalized 
by the Dirac form factor at zero momentum transfer, Zy = 1/-Fi(0). This renormalization 
is valid due to the good chiral properties of DWF. At zero momentum transfer, the result 
at rrif = 0.005 is smaller than the other masses which corresponds to the bending of 
discussed in Sec. IIV A[ Furthermore, the q 2 dependence of the results at the lightest quark 
mass is milder than the other masses. 

In the following we focus only on the momentum transfer dependence of the axialvec- 
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FIG. 20: Ratios of 2- and 3-point functions for axialvector current, A^(qs = 0,t),A^(q^ / 0, t) 
and A^(q,t) at m/ = 0.01. 
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FIG. 21: The axialvector form factor, Fa{q 2 ), renormalized by Zy = 1/Fi(0). 
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FIG. 22: The axialvector form factor, Fa{<1 2 ), normalized at q 2 = 0. The dashed line denotes a fit 
to experimental data. 

tor form factor: We normalize the form factor by its value at zero momentum transfer 
respectively for each quark mass. Figure [22] shows the results after these normalizations, 



FA_(q 2 ) / Fa(0) . For the heavier three masses, the results tend to decrease with quark mass 
while the dependence is opposite for the lightest mass. Similar to the vector-current form 
factors, the experimental axialvector form factor is also traditionally considered to be fitted 
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FIG. 23: Effective dipole mass M^ ff of FA(q 2 ) at mj = 0.005, 0.01 together with the experimental 
result Result of Dipole fit (solid line) with one standard deviation (dashed line) is also 

presented. 



well by the dipole form, 



1 



*U(o) 



(28) 

for the axialvector 



with the experimental data giving a best fit of Ma = 1.03(2) GeV [l^] 
dipole mass. The experimental fit is shown by the dashed line in Fig. [22j 

If the dipole form is valid in the entire q 2 region, we can extract the effective axial dipole 

mass, 

Mf = q2 = . (29) 



\ y/F A {0)/F A (?) - 1 

at each non-zero value of q 2 . Figure [231 shows that the effective dipole mass at m; = 0.01 is 
reasonably flat. This means that the form factor behaves as a dipole, as in the cases of the 
Dirac and Pauli form factors. We fit the form factor with the dipole form, and the fitted 
dipole mass is consistent with the effective one, as shown in Fig. [23] by the solid line with 
the one standard deviation (dashed lines). Figure [231 shows that the lightest quark mass 
data is also well explained by the dipole form, although the results do not approach the 
experimental value. 

The axial rms radius is determined from the dipole mass, 



(r A ) 



2\l/2 



12/M^ 



(30) 
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FIG. 24: Axial charge rms radius (r 2 ^) 1 ^ 2 determined from dipole fit. Dashed line represents a 
linear extrapolation of our data excluding the lightest point(striped circle). Square, up triangle, 

2] and quenched DWF and two- flavor and quenched 



diamond and left triangle denote two- flavor 



Wilson 



calculations, respectively. The star denotes the experimental result 



m f 0.005 0.01 0.02 0.03 m* s 

{r 2 A ) 1 ' 2 [fm] 0.366(36) 0.469(21) 0.423(25) 0.413(17) 0.493(33) 

g nNN (def.) 8.53(82) 10.38(94) 11.1(1.3) 12.0(1.1) 9.5(1.6) 

9nNN(GT) 11.84(45) 13.12(47) 12.66(78) 13.56(57) 12.79(79) 

g P 6.71(60) 8.45(71) 10.31(88) 11.93(93) 6.6(1.2) 

TABLE V: Axial charge rms radius (r 2 ^) 1 ^ 2 , nucleon-pion coupling g-nNN and induced pseudoscalar 
coupling gp. g-^NN is calculated with the definition eq. (|34ft and Goldberger-Treiman (GT) relation 
eq. (i35j) denoted as def. and GT in table, respectively. The linear fit results at mP hys = 135 MeV 
obtained without the lightest quark mass, are also presented. 

and is 0.666(14) fm in the experiment. The calculated axial rms radius from the fits is shown 
in Fig. [24] plotted as a function of the pion mass squared. The results are summarized in 
table|Vl While the result increases as the pion mass decreases, the lightest result significantly 
decreases. This pion mass dependence is similar to that observed in the axial charge in Fig. [71 
This, however, is not clear in F^q 2 ) renormalized by Zy = 1/^(0) obtained on our smaller 
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FIG. 25: Comparison of Fa(q 2 ) renormalized by Zy = 1/^(0) with larger and smaller volumes 
denoted by closed and open symbols, respectively, at each quark mass. Open symbols at q 2 = 
are slightly shifted to minus direction in x-axis. The dashed curve is a fit to experimental data. 

volume as shown in Fig. [25) the data at the lightest quark mass on the smaller volume shows 
a significant deviation from the larger volume result, but the statistical errors are too large 
to allow for a more quantitative comparison. 

Here, we similarly suspect this behavior of the larger volume to be caused by a large finite 
volume effect. A similar behavior is also seen in previous two-flavor results as presented in 
Fig. EH DWF [2] and Wilson 27] fermion calculations on a smaller volume (1.9 fm) 3 have 
similar pion mass dependences, but the radius begins to decrease at heavier pion mass. 
Once again, this behavior is quite similar to the case of the axial charge. Moreover, previous 
results obtained on large volumes jl|, Q do not exhibit such strong pion mass dependence, 
which is also shown in Fig. [2U Figure [26] shows the same results of the rms radii, but 
plotted as a function of m n L. The scaling of the rms radius with m n L is not as compelling 
as the axial charge case, but from the figure we estimate that m n L > 6 is required to obtain 
the axial charge rms radius without significant finite volume effects. Needless to say, other 
systematic errors, e.g., due to heavier quark mass than the physical one, should be removed 
to reproduce the experimental value. 

The lightest pion mass data is omitted in the following chiral extrapolation, because 
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FIG. 26: Same as Fig. [23] except the horizontal axis is the scaling variable m^L. Dashed lines 



denote the experimental result 



141 ] and its one standard deviation. 



we cannot rule out a large systematic error due to the finite volume of the simulations as 
discussed above. A linear fit to the heaviest three quark masses and extrapolation to the 
physical pion mass yields {r 2 ^) 1 ^ 2 = 0.49(3) fm. The fit result is presented in Fig. EU and 
reproduces 73% of the experimental value. 



2. Induced pseudoscalar form factor Fp(q 2 ) 

The induced pseudoscalar form factor, F P (q 2 ), is expected to have a pion pole, so its 
momentum-transfer dependence should be different from the other form factors. At the 
lightest quark mass this form factor is suspected to have a large finite volume effect, since it 
is obtained from the matrix element of the axialvector current together with the axialvector 
form factor, as discussed in the previous subsection. 

Figure [271 shows 2M^Fp{q 2 ) renormalized with Zy, plotted against the momentum trans- 
fer squared at each quark mass. We immediately notice that this form factor has a much 
larger q 2 dependence than the other form factors. In additi on, the results from all quark 



masses appear to be consistent with the experimental data 15[. Note that our statistical 
error is much smaller than the experiment. 

The induced pseudoscalar form factor is related to the axial vector form factor through the 
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FIG. 27: 2MjyFp(q 2 ) renormalized by Zy = 1/^(0) with experimental values 



so-called partially conserved axialvector current (PCAC) relation which is a manifestation 
of spontaneously broken chiral symmetry. In the traditional PCAC current algebra with 
pion-pole dominance (PPD), the PPD form, 

2M N F A (q 2 ) 



V) 



q 2 + m 2 



(31) 



is obtained at pa 0. The denominator on the right-hand side of this relation corresponds 
to the pion pole. We investigate the validity of this relation in our results through a quantity, 

(q 2 + ml)F P (q 2 ) 



(32) 



2M N F A (q 2 ) ' 

If the relation holds we obtain unity for this quantity at all q 2 . Figure [28] shows a PPD 
calculated using our lattice results for Fa and Fp. There is no significant q 2 dependence, 
and while the values are close to unity, they are systematically less than one. We fit these 



results by a constant for each quark mass, whose results are presented in 



15 



igJUand table ED 
651 ] within the 



16 



While all the fit results are consistent with the experimental data 
larger error of the experiments, they are about 10-20% smaller than the prediction of the 
PPD form. 

We should note that the quantity at the lightest quark mass looks similar to the others, 
but F A (q 2 ) at rrif = 0.005 is suspected to have large finite volume effect as discussed in 
the last subsection. This means that F P (q 2 ) at rrif = 0.005 is expected to suffer from a 
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FIG. 28: Induced pseudoscalar form factor Fp(q 2 ) normalized by nucleon mass, Fa(<1 2 ), and pion 



pole, along with experimental values 



0.005 



0.01 



0.02 



0.03 



a PPD 0.833(25) 0.837(29) 0.903(32) 0.873(31) 
[mPp^/m^ 2 1.044(39) 1.009(19) 0.940(24) 0.977(20) 

TABLE VI: a PPD and m^? le . 

similarly-sized effect at the same quark mass. Thus, it appears that the two large finite 
volume effects cancel in this ratio. 

We check the consistency of the pole mass in Fp(q 2 ) with the measured pion mass at 
each quark mass by observing that the pole mass is given by 

2 2a ppn M N F A (q 2 ^ 



K lc )' 



- q 



(33) 



F P (q 2 ) 

where we use the fact that a PPD ^ 1 in our data. Figure [291 shows that the ratio m F °^ c /m n 
is reasonably consistent with unity, and has no large q 2 dependence except for the lightest 
quark mass point, which has large statistical error. The values obtained from a constant fit 
are presented in table ED This consistency suggests that F P (q 2 ) does indeed have a pion 
pole structure, which is consistent with the PPD form, however a PPD 7^ 1 in our data. We 
confirmed that m F °^ and a PPD obtained from a monopole fit of IM^F^j Fp are reasonably 
consistent with the above results, but have larger errors. 
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FIG. 29: Ratio of the pole mass in Fp(q 2 ) and measured pion mass. 

The pion-nucleon coupling is related to the induced pseudoscalar form factor via the 
relation 



Q-kNN 



lim 

r.2 1. TT) 2 



{q 2 + ml)F P (q 2 
2F„ 



(34) 



where F n = 92.4 MeV. Combining the above relation with the PPD form, eq. (1311) . we obtain 
the Goldberger-Treiman (GT) relation 



gnNNF-x — M^gA- 



(35) 



In this relation we assume Fa{0) ~ Fa^—ttv^). As such it suffers from a small mismatch in 
momentum transfer. Nevertheless, if we substitute the experimental values for the quantities, 
we obtain g n NN — 12.9. 

Figure I3TJ1 shows two calculations for the nNN coupling, g n NN- one uses the definition 
of g-wNN and another the GT relation at each quark mass, plotted against the pion mass 
squared. In determining g n NN, we use the measured pion decay constant at each quark mass 



from Ref. 



we 



471 | . Results for g^NN from both methods are given in table |Vj From Fig 
observe that g n NN obtained from both methods displays only a mild m 2 dependence, with 
the exception of the lightest mass results which show a significant downward shift away from 
the trend set by the three heavier mass values. This of course is another manifestation of 
the large finite size effect observed in the axial charge (see Sec. IIV A|) . Hence, for the chiral 
extrapolation we simply employ a linear fit form and exclude the lightest mass point. We 
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FIG. 30: Two measurements of the ttNN coupling with Goldberger-Treiman relation and defini- 
tion of g n NN- The experimental value 66(] is indicated by the star. Dashed lines present linear 
extrapolations of our data without lightest point. 

obtain the results at the physical pion mass, g n NN — 9.5(1.6) from the definition eq. fl34l) . 
and g-xNN = 12.8(8) from the GT relation eq. (l35|) . The value obtained using the GT relation 
agrees with a recent estimation of the coupling g n NN = 13.3(9) obtained from forward 7riV 
scattering data [66j, and also with the previous result, g^NN — 11-8(3), from a quenched 



simulation performed using the Wilson action 



271 ] estimated by the GT relation. The 



result from the definition eq. .on the other hand, is consistent with a quenched DWF 



) tl) ' 

Rigorously speaking, the GT relation is not valid in our data, since our data do not satisfy 
the PPD form due to a PPD ~ 0.85. Thus, the difference between the two determinations of 
9-kNN can be explained by tt PPD . Further study of the GT relation is an important future 
work, since the relation should be satisfied in the chiral limit, and at zero momentum transfer. 

The induced pseudoscalar coupling for muon capture on the proton, gp = m fl Fp(q^) 
where (f c = 0.88m^, is defined with the muon mass and the induced pseudoscalar form 
factor Fp at the specific momentum transfer where the muon capture occurs, — > n+u^. 

Since Fp(q 2 ) has significantly large pion mass and momentum-transfer dependences due 
to the pion pole, we subtract this contribution before performing the momentum transfer 
and chiral extrapolations. To do this, we first define the quantity with pion-pole subtraction 
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tion. The experimental result [16|] is indicated by the star. The dashed line represents a linear 
extrapolation of our data excluding the lightest point. 



by 

F P (q 2 ) = (q 2 + ml)F P (q 2 ) (36) 

at each q 2 , and then extrapolate this to the required momentum transfer q 2 . The induced 
pseudoscalar coupling is estimated by the normalization factor with the physical pion mass 

Fp(q 2 c ) 



9 P 



m. 



q 2 + {m% 



phys\ 2 ' 



(37) 



where mP hys = 135 MeV, at each quark mass as shown in Fig. [3U The figure shows that 
the result is almost linear as a function of the pion mass squared and decreases toward the 
experimental result for the three heavier mass values. Again the lightest mass result is an 
exception caused by the finite volume effect in g^, as discussed in Sec. IIV Al The result at 
the physical pion mass, presented in table |V] is obtained from a linear fit to the heaviest 
three pion masses, and is consistent with the recent experiment [l6] and analysis 65]. This 
result also agrees with the previous quenched DWF result while it disagrees with a 



quenched Wilson determination 



271 ] 671 ] which is almost half of the experimental value. 
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V. CONCLUSIONS 



We have studied the isovector nucleon form factors with Nf = 2 + 1 flavors of dynamical 
quarks using the domain wall fermion action at a lattice cutoff of a -1 = 1.73 GeV. The form 
factors are calculated with four light quark masses, corresponding to a lightest pion mass, 
m w = 0.33 GeV, and with momentum transfers down to q 2 ~ 0.2 GeV 2 . 

We have found the axial charge decreases significantly at the lightest quark mass point 
on the larger volume while the effect sets in for heavier quark mass on the smaller volume. 
By comparing our results with those using different volumes, numbers of flavors, and lattice 
fermions as a function of the single variable m^L, we conclude that this downward trend is 
caused by the finite volume used in our calculation. The fact that such an effect is absent in 
quenched and partially-quenched mixed-action studies on large volumes may be explained 
by the presence of unphysical logarithms. 

We have fit the data to several forms, including finite volume effects, and obtain g& = 
1.19(6)(4), where the first and second errors are statistical and systematic, respectively, 
which is 7% smaller than the experimental value. In our estimation, a spatial volume of 
V > (3.5 fm) 3 is required to keep the finite volume effect at, or below, one percent at 
m n = 0.33 GeV. Hence lattice calculations should continue to push down the quark mass 
and increase the volume with m n L > 6. Detailed analyses of the quark mass and finite 
volume dependence is desirable to understand the systematic deviation from experiment. 

Our lattice results for the form factors of the vector current are well fit by the standard 
dipole form. We have evaluated the root-mean squared radii, and the difference of the 
anomalous magnetic moment between the proton and neutron from the dipole fits. The 
radii and the anomalous moment are well explained by a linear function of the pion mass 
squared. In the radii we have not observed divergent quark mass behaviors predicted by 
HBChPT. Besides the divergent behavior, the pion mass dependences for the observables 
are quite consistent with other lattice QCD calculations including the recent results of LHP. 
Due to the linear behavior, we have concluded that the form factors of the vector current 
is less sensitive to the finite volume effect in contrast to the axial charge. Although both 
Dirac and Pauli rms radii approach to the experimental values as the pion mass decreases, 
the values extrapolated by the linear form at the physical pion mass underestimate the 
experiments by about 25%. Future work will involve simulating at lighter quark masses to 
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search for the nonanalytic behavior predicted by HBChPT. 

The axialvector form factor is also well described by the dipole form, even at the lightest 
quark mass, where the axial charge, -Fa(O), is suspected to have a large finite volume effect. 
The axial charge radius, obtained from the dipole fit, has a downward tendency as a function 
of the pion mass squared, which drives the radius away from the experimental value. We 
have considered this dependence to be caused by the finite volume of our simulation, as 
in the case of the axial charge. We observe that our results seem to scale do 
previous calculations using several volumes. We have concluded that the form factors of the 
axialvector current are more sensitive to the finite volume than those of the vector current 
from the observations of the finite volume effects. 

We have checked the pion-pole structure in the induced pseudoscalar form factor with 
our simulations. We have found that the pion-pole dominance form describes our data well, 
with the exception that a PPD < 1. Taking into account that a PPD 7^ 1, the pole mass of the 
induced pseudoscalar form factor reasonably agrees with the measured pion mass. 

For a precision test of QCD from nucleon matrix elements, we have identified several 
problems that first need to be overcome, such as finite volume systematic errors in the axial 
charge and the form factors of the axialvector currents, and the underestimation of the radii 
of the form factors of the vector current. Further lighter quark mass and larger volume 
calculations are essential to solve the problems, and such simulations are underway. Besides 
the comparisons with the experimental values, it is also important future work to study why 
a PPD deviates from unity. 
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m,f 


q 2 [GeV 2 l 


FAq 2 ) 


F 2 (q 2 ) 


F A {q 2 ) 


2M N F P (q 2 ) 


0.005 


0.0 


1.0000(11) 


N/A 


1.073(39) 


N/A 




0.198 


0.785(19) 


2.20(13) 


0.959(34) 


13.32(91) 




0.383 


0.622(22) 


1.716(97) 0.892(39) 


7.42(51) 




0.557 


0.505(28) 




U. i O^t^l j 


5.27(42) 




0.723 


0.516(53) 


l.OUl ID ) 


n 7Q r >((\K\ 
u. * yz^oo j 


4.71(53) 


0.01 


0.0 


1.0000(10) 


IV / A 

IN / 1\ 


l.lou^oo j 


N/A 




0.199 


0.787(17) 






13.00(90) 




0.385 


0.641(22) 






7.48(44) 




0.562 


0.524(31) 


1 1A(~\ 1 \ 


n 71 Q^CA 


5.01(42) 




0.731 


0.506(49) 




n 7m 


3.95(46) 


0.02 


0.0 


1.0000(15) 


N/A 


1.174(37) 


N/A 




0.200 


0.805(20) 


2.40(15) 


1.005(33) 


12.3(1.2) 




0.390 


0.686(32) 


2.08(13) 


0.890(38) 


9.04(73) 




0.573 


0.599(49) 


1.80(13) 


0.839(48) 


7.06(70) 




0.748 


0.443(37) 


1.31(14) 


0.668(47) 


3.84(59) 


0.03 


0.0 


1.0000(11) 


N/A 


1.196(30) 


N/A 




0.201 


0.8302(99) 


2.79(12) 


1.038(28) 


12.74(88) 




0.394 


0.700(15) 


2.302(92) 0.912(32) 


9.68(62) 




0.580 


0.595(22) 


2.00(12) 


0.838(40) 


6.81(61) 




0.760 


0.500(31) 


1.54(15) 


0.704(50) 


4.55(73) 



TABLE VII: Form factors of vector, axialvector currents on (2.7 fm) 3 . All form factors are renor- 
malized. 
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m f q 2 [GeV 2 ] F.tf) F 2 (q 2 ) F A (q 2 ) 2M N F P (q 2 ) 



0.01 


0.0 


1.000(27) 


N/A 


1.066(72) 


N/A 




0.430 


0.621(52) 


1.42(23) 


0.580(93) 


3.2(1.4) 




0.812 


0.46(14) 


0.89(50) 


0.380(70) 


1.83(68) 


0.02 


0.0 


1.000(14) 


N/A 


1.115(58) 


N/A 




0.437 


0.705(38) 


2.03(17) 


0.749(49) 


6.16(97) 




0.833 


0.501(36) 


1.14(14) 


0.504(47) 


2.70(38) 


0.03 


0.0 


1.0000(6) 


N/A 


1.149(32) 


N/A 




0.441 


0.686(18) 


1.91(10) 


0.787(29) 


6.56(70) 




0.848 


0.522(23) 1.258(70) 0.574(30) 


3.48(30) 



TABLE VIII: Form factors of vector, axialvector currents on (1.8 fm) 3 . All form factors are 
renormalized. 
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